Multi-band Superconductivity in the Chevrel Phases SnMoeSg and PbMoeSs 
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Sub-Kelvin scanning tunnelling spectroscopy in the Chevrel Phases SnMoeSs and PbMoeSs re- 
veals two distinct superconducting gaps with Ai = 3 meV, A2 ~ 1.0 meV and Ai = 3.1 meV, 
A2 ~ 1.4 meV respectively. The gap distribution is strongly anisotropic, with A2 predominantly 
seen when scanning across unit-cell steps on the (001) sample surface. The spectra are well-fitted 
by an anisotropic two-band BCS s-wave gap function. Our spectroscopic data are confirmed by 
electronic heat capacity measurements which also provide evidence for a twin-gap scenario. 



PACS numbers: 

Among the vast zoo of poorly-understood supercon- 
ductors, Chevrel Phases (CP) stand out for their high 
upper critical fields Hc2, niany of which exceed the 
Pauli limit These materials were first synthesised 
in 1971 and enjoyed a wealth of attention in the early 
1980s. Unfortunately, the discovery of the cuprate super- 
conductors largely swept CP under the laboratory carpet, 
despite a lack of detailed understanding of their large 
Hc2 values. A multi-band scenario (incorporating strong- 
coupling effects and enhanced spin-orbit scattering) was 
suggested as a possible explanation "3" , but until now this 
hypothesis has remained experimentally unexplored. 

Multi-band superconductivity was first proposed 50 
years ago as a potential avenue for increasing critical 
temperatures g]. Interband scattering between non- 
degenerate bands at the Fermi level Ep enables supercon- 
ductivity to be induced in bands which may not directly 
participate in the pairing mechanism, thus increasing the 
effective density of states (DoS) and hence the transi- 
tion temperature T^. However, with the exception of 
some transition metal calorimetric data [S] and tunnelling 
in doped SrTiOa p|, multi-band superconductivity re- 
mained a largely theoretical concept until the discovery of 
MgB2 in 2001 revived interest in the field [7]. In this ma- 
terial, superconductivity in the quasi-2D tr-band induces 
coherence in the quasi-3D 7r-band with an unexpectedly 
high Tc of 39 K. The two gaps have been imaged by a va- 
riety of techniques, including local spectroscopic [H] and 
bulk thermodynamic approaches Recently, evidence 
has been found for multi-band superconductivity in boro- 
carbides [inj, sesquicarbides [TT], skutterudites [T^] and. 



perhaps most interestingly, pnictides [13 . CP and pnic- 
tides share similar anomalously large values of Hc2 and do 
not follow standard Werthamer-Helfand-Hohenberg the- 
ory. However, in contrast with the pnictides, the MogXg 
{X — S, Se) Chevrel cluster does not exhibit any intrin- 
sic magnetism or competing order. This greatly simpli- 
fies the analysis and interpretation of its low-temperature 
properties, particularly any multi-band effects. Band- 
structure calculations have indicated the presence of two 
Mo d bands at Ep in CP [Mj: in this Letter we present 
local spectroscopic evidence for two distinct supercon- 
ducting gaps in SnMogSg and PbMogSg. These data are 
supported by specific heat measurements displaying clear 
signatures of a second gap. 

We have chosen to focus on SnMogSg and PbMogSg 
since these two materials have the highest values for Tc 
and Hc2 within the CP family: 14.2 K, ~ 40 T and 
14.9 K, > 80 T respectively [15]. Single crystals of each 
compound with typical volume 1 mm^ were grown at 
1600 °C by a chemical flux transport method using sealed 
molybdenum crucibles. Their high purity was confirmed 
by AC susceptibihty (ACS) yielding AT^ = 0.1 K for 
SnMoeSg and 0.3 K for PbMogSg. Local spectroscopy 
(STS) was performed on room-temperature-cleaved sam- 
ples with a home-built helium-3 scanning tunnelling mi- 
croscope in high- vacuum (< 10~^ mbar), using a lock-in 
amplifier technique. Heat capacity measurements were 
carried out at the Grenoble High Magnetic Field Labo- 
ratory with a high-resolution microcalorimeter using the 
"long relaxation" technique |16j and in Geneva using a 
Quantum Design"^^ PPMS. 
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The first hint of a two-band order parameter arises 
from fast spectroscopic traces over several tens of 
nanometers in the (001) plane of each material (Fig. [ij. 
The corresponding topography in SnMogSg shows atom- 
ically flat terraces separated by steps of size 12 ± 1 A, 
which compares favourably with twice the rhombohedral 
unit cell parameter 6.5 A. Spectra taken on the terraces 
are homogeneous, with a gap of 2.95 meV and a marked 
lack of any quasiparticle excitations within the gap. In 
contrast, spectra taken on the steps between terraces dis- 
play additional peaks at lower energy, suggestive of a 
second gap. We interpret this as a local modification of 
the tunnelling matrix element, enabling us to preferen- 
tially probe another portion of the Fermi surface with 
different atomic orbital characters. Cleaved surfaces of 
PbMogSs are of rather lower quality with an RMS rough- 
ness of ~ 1.5 A and broad poorly-resolved unit cell-sized 
steps (Fig. [TJb)(iii)). However, the average spectrum 
(Fig. [ljb)(ii)) displays a kink at ~ ± 1.4 meV (high- 
lighted by arrows) and a V-shaped dispersion around Ep. 
This confirms the presence of states within the large gap. 

Such a dramatic spectral variation as a function of the 
local topography has not previously been observed in any 
other superconductor. It may therefore be natural to sug- 
gest that the isolated appearance of these multi-gap sig- 
natures at unit cell steps could be due to a surface bound 
state or defect. However, a localised state would not dis- 
play the particle-hole symmetry of the peaks we observe. 
We have imaged a large number of separate unit cell steps 
and a second gap is consistently observed upon scanning 
across them. Another explanation for the double-gap 
behaviour could be the proximity effect inducing weak 
superconductivity in a metallic surface layer il7J. How- 
ever, the small gap induced would vary strongly with the 
thickness of the surface metallic layer. Apart from the 
fact that measurements are performed on freshly-cleaved 
samples, thus rendering any surface layer deposition im- 
plausible, a layer of metallic impurities would not be ex- 
pected to have a uniform thickness. This would cause 
substantial variation in the size of the induced gap and 
an extremely high zero-bias conductance (ZBC), both of 
which are incompatible with our data. 

In Figure [2] we display a range of spectra with fits us- 
ing a multi-band model. The Bardeen-Cooper-Schrieffer 
(BCS) quasiparticle density of states for an anisotropic 
s-wave n-band superconductor may be written as 



Rc 



(w + iTj)sign{uj)d0 



(1) 



where Nj is the contribution of band j to the DoS at 
Ep: Fj the scattering rate due to lifetime effects, Aj 
the magnitude of the gap within band j and Fj(9) — 
flj + (1 — aj) cos 9 measures the anisotropy of the corre- 
sponding gap with 0.5 < aj < 1. We include the temper- 
ature and the experimental smearing (0.3 meV) before 
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FIG. 1: (a) Zero-field 35nm trace on SnMoeSg taken at T 
= 0.4K, junction resistance Rt = 0.03 GO., (i) Topography 
showing double unit-cell steps; (ii) spectroscopic trace; (iii,iv) 
individual spectra taken on a flat terrace (1) and above a unit- 
cell step (2). (b) Zero-field 40nm trace on PbMoeSg taken at T 
— 0.5K, Rt = 0.015 GQ,. (i) Spectroscopic trace; (ii) average 
spectrum from entire trace; (iii) topographic variation. All 
data are raw and unaveraged. 



performing least-squares fits to our data with Nj , Aj , Fj 
and Vj as free parameters. Note that the spectral back- 
grounds between ± 5-10 meV are rather poorly-fitted, 
which is indicative of strong coupling to a low-energy 
phonon. However, a full Eliashberg analysis of the spec- 
tra is beyond the scope of this Letter. 

Atomically flat surfaces in SnMogSg produce homoge- 
neous spectra (Fig[2f^a)) which may be fitted using only 
a single band (i.e. n=l in ([!])). There is a slight deterio- 
ration in the fit quality at low energy, which is attributed 
to a very small contribution from the second band. In 
contrast. Fig. [2jb) shows the average of around 50 spec- 
tra acquired above a unit cell step. There is clearly a 
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FIG. 2: (a-c) SnMoeSs spectra and fits taken at 0.4 K, Rt 
= 0.03 GO., (d-e) PbMoeSg spectra and fits taken at 0.5 K, 
Rt = 0.015 Gn. (f) PbMoeSg spectrum taken at 1.9 K, Rt 
= 0.025 GO. from [E]- See text for details and table |l] for fit 
parameters. 



significant contribution from the smaller gap, necessi- 
tating a two-band fit. Similar fits are carried out on 
spectra from a flat zone and a broad step in PbMoeSs 
and the parameters obtained listed in Table |l] We find 
2Ai/fcBTc ~ 5 in each compound, but A2 is 30-40 % 
larger in PbMoeSg than SnMogSg. In both materials the 
gap anisotropies are similar: a small anisotropy in Hc2 
(e^ = 0.67) has been observed in PbMoeSs [3], but with 
the present data we are unable to judge whether this is 
due to the anisotropy in Ai or A2. We believe it unwise 
to draw quantitative conclusions on the symmetry of A2, 
since our experiment has a finite resolution imposed by 
a 0.3 meV broadening from the lock-in. However, any 
interband scattering will preclude a pure d-wave order 
parameter in A2 due to the dominant isotropic s-wave 
component in Ai. 

Previous STS experiments on PbMogSs provided ev- 
idence for low-energy excitations within the supercon- 
ducting gap, but lacked sufficient resolution to distin- 
guish two separate gaps. This is due to three factors: 
sample age, temperature and environment. In |18j . mea- 
surements were performed on old crystals at 1.9 K in an 
exchange gas, compared with freshly-grown samples at 
0.4-0.5 K and high vacuum in the present work. The in- 
creased thermal broadening at 1.9 K blurs the two gaps, 
though this should not be sufficient to render the smaller 
gap invisible. The major factor here is a deterioration 
in the sample surface due to the exchange gas environ- 
ment. It is well-known that in a two-band superconduc- 
tor, interband scattering due to impurities mixes the two 
gaps and reduces T^, resulting in an effective single-band 
anisotropic superconductor in the dirty limit. This was 
first predicted for MgB2 [121 [2D] and later observed in 



TABLE I: Superconducting gap parameters and relative DoS 
contributions from tunnelling (STS), heat capacity (HC) and 
AC susceptibility (ACS) data. Both gaps Ai,2 are measured 
in meV. ri,2 < 0.2 meV for all STS fits. 





SnMogSs 


PbMogSg 


technique 


bulk 


bulk 




ACS 


14.2 ± 0.05 K 


14.9 ± 0.15 K 




HC 


42 ± 1 T 


86 ± 5 T 


7 


HC 


6.4 ± 0.1 m 


Jgat~^K 


6.7 ± 0.1 mJgaf^K"^ 


Al 


HC 


3.06 


± 0.1 


3.15 ± 0.1 


A2 


HC 


0.86 ± 0.1 


1.41 ± 0.1 


Ni 


HC 


96 ± 2% 


90 ± 2% 


N2 


HC 


4 ± 2% 


10 ± 2% 




Terrace 


Stop 


Terrace 


Step 


Al 


STS 


2.92 ± 0.1 


2.95 ± 0.1 


3.14 ± 0.15 


3.06 ± 0.15 


A2 


STS 




1.05 ± 0.2 


1.42 ± 0.2 


1.36 ± 0.2 


ai 


STS 


0.85 ± 0.02 


0.87 ± 0.02 


0.85 ± 0.02 


0.89 ± 0.02 


0.2 


STS 




0.91 ± 0.1 


0.92 ± 0.1 


0.75 ± 0.1 


Ni 


STS 




62 ± 4% 


90 ± 4% 


66 ± 4% 


N2 


STS 




38 ± 4% 


10 ± 4% 


34 ± 4% 



irradiated samples [21| . However, due to extremely weak 
scattering between a and tt bands, the single band limit 
is never reached in MgB2. This may not be the case 
for CP: Fig. |2jc) displays a SnMogSs spectrum from a 
terrace after 3 months of measurements comprising nu- 
merous thermal and magnetic cycles. It is qualitatively 
similar to the results of [18i (shown in Fig. [2]jf)), provid- 
ing good evidence for low-energy states within the large 
gap, but does not display a distinct smaller gap. This 
is consistent with the presence of strong interband sur- 
face scattering. The ZBC is also rather high in both (c) 
and (f), which we attribute to a decrease in the super- 
fluid density due to enhanced pair-breaking from inelastic 
scattering. 

Upon increasing the temperature the large gap is grad- 
ually reduced and closes at the bulk determined by 
ACS. No pseudogap is visible above T^, conflrming that 
superconductivity arises from a metallic ground state and 
hence justifying the use of a BCS model to fit the spec- 
tra. In Fig. [3] we have plotted the variation of the large 
gap Al with temperature for each compound, with the 
theoretical BCS weak-coupling s-wave curve for compar- 
ison. A small kink is visible within each curve (shaded 
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FIG. 3: Variation of Ai(T) in (a) SnMoeSg and (b) PbMoeSg. 
The gap was determined by fitting with a BCS single-band 
anisotropic s-wave model for simplicity and all spectra were 
acquired on a flat terrace. 
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areas). Similar features have been observed for the tt- 
band (smaller) superconducting gap in MgB2 and origi- 
nate from interband scattering "stretching" the effective 
Tc of the weakly-coupled tt band to the bulk T^. We 
suggest that in CP a small contribution from a weakly- 
coupled band "stretches" the intrinsic Tc of a strongly- 
coupled band (which provides the majority of the DoS) 
to the measured bulk Tc. The position of the kink at 
higher energy and lower temperature in PbMogSg com- 
pared with SnMogSg is consistent with our observation 
that band 2 is more strongly-coupled in FbMogSs. We 
hypothesise that this may be the key to PbMogSs hav- 
ing a significantly higher Hc2 than SnMogSg, although 
further experiments will be required for confirmation. 

It is instructive to complement our STS measurements 
with bulk thermodynamic (HC) data, in order to con- 
clusively rule out any spurious surface effects being re- 
sponsible for A2 . Figure |4] (a) and (b) display the elec- 
tronic heat capacity Cdoc in SnMogSg and PhMogSg: this 
is measured by subtracting the HC in an applied field 
H = 28 T from the zero-field data Cqt- To eliminate 
the effect of fluctuations above Tc in high field, we limit 
our data to T > 1.15 rc(28T), where rc(28T) = 4.4 K 
and 9.3 K in SnMogSg and PhMogSs respectively. In 
CP there is a large contribution to the lattice HC Ciatt 
from a low-energy Einstein phonon due to vibrations of 
the cation between the MogXg clusters: it is therefore 
not possible to calculate Ciatt and hence extract Cdcc 
using the conventional S"^]^ i?2j+iT^^+^ model. How- 
ever, the Sommerfeld constant 7 may still be calculated 
using entropy considerations (see Table |l]). We find that 
■'j/Tc = 0.45 mJgat^-^K"'^ in each compound, suggesting 
that Tc scales with the DoS at Ep. Using a two-band 
a-model |22j we have performed fits to Coioc and sum- 
marise our results in Table |l] For both gaps, the values 
of 2A/kBTc from our STS data agree perfectly with those 
from HC experiments. While it is not possible to quan- 
titatively compare our STS-measured Nj (which also de- 
pends on the tunnelling matrix element) with the bulk 
Nj, the trends observed by each technique (iVi>iV2) are 
qualitatively in agreement. 

The final signature of 2-band superconductivity is pro- 
vided by the low-temperature variation of Ccicc{H) in 
each material. In a single-band BCS s-wave supercon- 
ductor ^{H) should be linear. However, at T = 0.35 K 
we observe bends in Ceiec{H) at = 2.8 ± 0.2 T and 
3.4 ± 1 T in SnMogSg and PbMoeSg respectively, remi- 
niscent of the low-field behaviour of MgB2. Extrapolat- 
ing the high-field linear fits to the normal-state value for 
7, we obtain 7Jc2 = 42 ± 1 T and 86 ± 5 T (although it 
should be noted that these may be slight overestimates 
due to vortex overlap effects at high field). We assume 
that Hx corresponds to the crossover between filling A2 
in band 2 followed by Ai in band 1 and hence estimate 
A^i=93 ± 0.5 %, N2=7 ± 0.5 % and iVi=96 ± 1 %, 
iV2=4 ± 1 % for SnMogSg and PhMogSg. These figures 
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FIG. 4: (a),(b) Cdcc/'yT with two-band a-model fits [22]. 
Insets: C/T at T and 28 T. (c),(d) Field-dependent con- 
tributions to Ccioc at T=0.35 K with linear fits above and 
below the crossover field (see text). The crossover in 
FbMogSg is rather broad compared to SnMoeSg; this is due 
to the lower homogeneity in the PbMogSg crystal as seen by 
increased transition widths in HC and ACS data. 



are in good agreement with those in Table [T] 

Together, our spectroscopic and thermodynamic data 
provide compelling evidence for a multi-band order pa- 
rameter in CP superconductors. In both SnMogSg and 
PbMogSg, a strongly-coupled quasi-isotropic band (con- 
tributing the majority of the DoS at Ep) coexists with a 
highly anisotropic weakly-coupled minority band. Look- 
ing ahead, we postulate that understanding and manipu- 
lating the interplay between two or more such bands may 
hold the secret to realising high values for Hc2 in future 
superconducting materials. 
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